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Abstract 

We study the metal-insulator transition (MIT) in weakly coupled disordered 
planes on the basis of a Non-Linear Sigma Model (NLerM). Using two differ- 
ent methods, a renormalization group (RG) approach and an auxiliary field 
method, we calculate the crossover length between a 2D regime at small length 
scales and a 3D regime at larger length scales. The 3D regime is described by 
an anisotropic 3D NLcM with renormalized coupling constants. We obtain 
the critical value of the single particle interplane hopping which separates the 
metallic and insulating phases. We also show that a strong parallel magnetic 
field favors the localized phase and derive the phase diagram. 

I. INTRODUCTION 

A quasi-two-dimensional (quasi-2D) weakly disordered conductor (i.e. weakly coupled 
planes system) exhibits a MIT for a critical value of the single particle interplane hopping t±. 
The existence of this quantum phase transition results from the localization of the electronic 
states in a 2D system by arbitrary weak disorder, while a 3D system remains metallic below 
a critical value of the disorder. 

Although the existence of this MIT is well established, its properties are far from being 

(c) 

understood. In particular, the critical value t ± of the interplane coupling remains contro- 
versial. The self-consistent diagrammatic theory of the Anderson localization! predicts an 
exponentially small critical coupling! t± greement with simple 

scaling arguments! and other analytical arguments.! (Here r is the elastic scattering time, 
kp the 2D Fermi wave vector and I the intraplane mean free path, kpl ^> 1 for weak 
disorder.) A recent numerical analysis predicts a completely different result, t± ~ \ j\fr, 
which is supported by analytical arguments.! It has also been claimed, on the basis of dia- 
grammatic perturbative calculations, that the MIT depends on the propagating direction, 
in contradiction with the scaling theory of localization.^! (Note that for a system of weakly 
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coupled chains, it is well established, both from numerical and analytical approaches, that 
t ( l ] ~ l/r.0) 

The aim of this paper is to reconsider the MIT in quasi-2D conductors on the basis of 
a NLaM.i In the next section, we use a renormalization group (RG) approach to show how 
the system crosses over from a 2D behavior at small length scales to a 3D behavior at larger 
length scales. The crossover length L x is explicitly calculated. We show that the 3D regime 
is described by an effective NLaM with renormalized coupling constants. We obtain the 
critical value of the interplane coupling and the anisotropy of the correlation (localization) 
lengths in the metallic (insulating) phase. Next, we show that a parallel magnetic field 
tends to decouple the planes and thus can induce a MIT. The results of the RG approach 
are recovered in section [TT1] by means of an auxiliary field method. We emphasize that the 
latter is a general approach to study phase transitions in weakly coupled systems. 

II. RENORMALIZATION GROUP APPROACH 

We consider spinless electrons propagating in a quasi-2D system with the dispersion law 

e k = kj/2m - 2t± cos(k ± d) , (2.1) 

where kii and k± are the longitudinal (i.e. parallel to the planes) and transverse components 
of k, respectively, m is the effective mass in the planes, d the interplane spacing and t± the 
transfer integral in the transverse (z) direction. The effect of disorder is taken into account 
by adding a random potential with zero mean and gaussian probability distribution: 

(^(r)ty(rO) = (27riV 2 (0)r)- 1 5 M ,5(r - r') , (2.2) 

where AT 2 (0) = m/27t is the 2D density of states at the Fermi level and r the elastic scattering 
time, r is the coordinate in the plane and the integers I, /' label the different planes. We 
note kp the 2D Fermi wave vector and vp = kpjm the 2D Fermi velocity. 



A. A simple argument 

We first determine the critical value by means of a simple argument whose validity will 
be confirmed in the next sections. Consider an electron in a given plane. If the coupling t± is 
sufficiently weak, the electron will first diffuse in the plane and then hop to the neighboring 
plane after a time r x . The corresponding diffusion length L x is determined by the equations 

L 2 x = D(L x )r x , 

d 2 = D ± t x = 2t\d 2 TT x . (2.3) 

The length dependence of the coefficient D(L) results from quantum corrections to the 
semiclassical 2D diffusion coefficient D = v 2 F r/2. We have assumed that the transverse 
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diffusion is correctly described by the semi classical diffusion coefficient D±. As shown in 
the next sections, this is a consequence of the vanishing scaling dimension of the field in the 
NLcrM approach. Eqs. ( |2.3| ) give t x ~ l/t 2 ± T and L 2 ~ D(L x )/tj_r. The critical value t± is 
obtained from the condition L x ~ £ 2 £> or g(L x ) ~ 1, where is the 2D localization length 
and g(L) the 2D (dimensionless) conductance. These two conditions are equivalent since 
g(L) = N 2 {0)D(L) ~ 1 for L ~ &d (see Eq. ( gT3|) below). We thus obtain 

tj c) ~ , 1 !=i, (2.4) 

>$ d N 2 (0)t £WMt 



where / = v^r is the elastic mean free path. In the weak disorder limit (kpl 3> 1), C.2D ~ 
le akpl (a ~ 1) so that t± is exponentially small with respect to 1/r. Apart from the 
factor 1 / y/kpl, Eq. (|2.4|) agrees with the result of the self consistent diagrammatic theory of 
Anderson localizations and with estimates based on the weak localization correction to the 
Drude-Boltzmann conductivity.! 



B. NLcrM for weakly coupled planes 



The procedure to derive the NLcrM describing electrons in a random potential is well 
established! and we only quote the final result in the quasi-2D case (see however section [11 D| ) . 
Averaging over disorder by introducing N replica of the system and using an imaginary time 
formalism, the effective action of the model is Sa = S 2 d + S± with 

S 2 d[Q) = lN 2 (0)DY: J dVTrfViiQKr)] 2 - |iV 2 (0) £ J d 2 rTr[^(i-)] , 
Sa.[Q] = -jN 2 (0)tlr ^ / rf 2 rTr[g,(r)gr(r)] , (2.5) 

where Qi(r) is a matrix field with elements ijQ?n m {r). a, (3 = 1...N are replica indices, n,m 
refer to fermionic Matsubara frequencies uu n ,u m , and i,j = 1,2 describe the particle-hole 
(i = j) and particle-particle (i 7^ j) channels. Q is the matrix ift^nm = ^i,j^n,mS a ,/3^n- in 
( p.5| ), Tr denotes the trace over all discrete indices. The field Q satisfies the constraints 
Qf(r) = 1 (with 1 the unit matrix), TrQ ; (r) = and Q + = C T Q T C = Q where C is the 
charge conjugaison operator.! A ~ l/l is an ultraviolet cut-off for the fluctuations of the 
field Q in the planes. The Josephson-like interplane coupling S± is obtained retaining the 
lowest order contribution in t± ((1,1') denotes nearest neighbors). 

If t^r ^> DA 2 ~ 1/r, the fluctuations in the transverse direction are weak. In this 
case, Tr[QiQi/} can be approximated by — (d 2 /2)Tr[V Z Q} 2 and we obtain a NLcrM with 
anisotropic diffusion coefficients D and D±. For t±r ^> 1, an electron in a given plane hops 
to the neighboring plane before being scattered. There is therefore no 2D diffusive motion 
in that limit so that the quasi-2D aspect does not play an essential role. The anisotropy can 
be simply eliminated by an appropriate length rescaling in the longitudinal and transverse 
directions. 00 We consider in the following only the limit t±r 1 where we expect a 2D/3D 
dimensional crossover at a characteristic length L x . 
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C. RG approach 



We analyze the action ( |2.5| ) within a RG approach following Ref. [12|. The condition 
t_i_r <C 1 ensures that the transverse coupling is weak (DA 2 ^> t\r). The initial stage of 
the renormalization will therefore be essentially 2D, apart from small corrections due to the 
interplane coupling. If we neglect the latter, we then obtain the renormalized action 

Sa>[Q] = ^iV 2 (0) J D(A')E / d 2 rTr[V||Q z (r)] 2 - |iV 2 (0) £ / d 2 rTr[^(r)] 

-jN 2 (0)tlr £ / d 2 rTr[Q,(r)Q,(r)] , (2.6) 

where A' < A is the reduced cut-off after renormalization and D(A') the renormalized 
value of the 2D longitudinal diffusion coefficient. Since the scaling dimension of the field Q 
vanishes (i.e. there is no rescaling of the field) ,@ there is no renormalization of the interplane 
coupling.li3 The 2D/3D dimensional crossover occurs when the transverse and longitudinal 
couplings become of the same order: 

l -D x hl~tlT, (2.7) 

where D x = D(A X ) is the longitudinal diffusion coefficient at the crossover. In the 3D regime 
(A' < A x ), it is appropriate to take the continuum limit in the transverse direction!!! Using 

T±[QiQi>] = -\^\Qi ~ Qv? + const - -^Tr[V z Q,] 2 (2.8) 
for I and /' nearest neighbors, and dJ2i I dz, we obtain 

SaM = lN 3 (0)D x J d 3 r[Tr[V||Q(r)] 2 + (dA x ) 2 Tr[V 2 Q(r)] 2 ] - ^iV 3 (0) J d 3 rTr[fiQ(r)] , 

(2.9) 

where AT 3 (0) = AT 2 (0)/(i is the 3D density of states at the Fermi level. The cut-offs are A x 
and 1/d in the longitudinal and transverse directions, respectively. Note that r is now a 
3D coordinate. The 3D regime is thus described by an anisotropic NLcrM. However, the 
anisotropy is the same for the diffusion coefficients and the cut-offs and can therefore easily 
be suppressed by an appropriate rescaling of the lengths: rj, = ry/si, z' = z/s 2 , with 

sfD x = sfD x (A x d) 2 , 

s 2 s 2 = l. (2.10) 



The last equation ensures that / <fV = / d 3 r and lets invariant the last term of ( j2.9| ). From 
( P-lOj ), we obtain si = (l/A^d) 1 / 3 , s 2 = (A x d) 2 ^ 3 , and the new effective action 

SaM = In 3 (0)D J d 3 rTr[VQ(r)} 2 - |iV 3 (0) J d 3 rTr[fiQ(r)] , (2.11) 
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where D = D x (A x d) 2/3 . The cut-off is now isotropic: SiA x ~ s 2 /d ~ A = (A 2 x /d) 1/3 . The 
dimensionless coupling constant of the NLcxM ( |2.11 ' 



is 



4 4 4 

A = 7r^ 3 (0)D A= 7rA^0)^ = ^/ (2 - 12) 

The MIT occurs for A = A c = 0(1), i.e. when the 2D conductance g x = g c = 0(1). Using 
9x — 9c ~ 1 for A x ~ we recover the result (|2.4| ) for the critical value of the interplane 
coupling. 

Further information about the metallic and insulating phases can be obtained from the In- 
dependence of the 2D dimensionless conductance. The self-consistent diagrammatic theory 
of Anderson localization givesS 



g(L) = -ln 1 + ^ 1 + - e~ L ^ , (2.13) 




where L should be identified with A" 1 . Eq. Q2.13Q gives g(L) ~ (1/tt 2 ) ln(£ 2D /L) for L < £ 2 d 



(in agreement with perturbative RG calculation^) and g(L) ~ (l/27r 2 )(^2D/-^)e^ L ^ 2D for 
L ^> ^2£>- The crossover length L x obtained from (|2.7| ) and ( |2.13| ) is shown in Fig. [I]. Since 
the crossover length is not precisely defined, we have multiplied L x by a constant in order to 
have L x ~ vp/t± deep in the metallic phase (t± 3> t±) (this will allow a detailed comparison 
with the results of Sec. PTT|) . 

In the metallic phase, t± > t^, we have 

^~2iV 2 (0)tir^ D e- 2 ^ ^1, 

~ £ 2D e"^ < 6d • (2.14) 



Since the renormalization in the 3D regime does not change the anisotropy, we obtain from 
( p.9| ) the anisotropy of the correlation lengths (see Fig. |2|): 

L\\ <r ± (A x d) 2 ' 1 • ' 

Deep in the metallic phase (t± ^> t^), we have L±/L\i ~ (t||/tj_) 2 (where t\\ ~ ■uf&f ~ vp/d 
is the transfer integral in the planes) while close to the transition (t± ^ t^) we have L±/L\\ ~ 
i&D/d) 2 ~ (^||A±) 2 / (&fO- Thus, as we move towards the transition, the anisotropy deviates 
from the result L±/Lu ~ (£||/tj_) 2 predicted by the 3D anisotropic NLctM.0 
In the insulating phase, t± < , we have 



1 , 

g x - ^20^ 



2Af 2 (0)tir 



-e 



gx 1 
2iV 2 (0)t^T (.2D < 1 



^«3? iS6D - (2 - 16) 



The anisotropy of the localization lengths is given by (see Fig. 
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Close to the transition, A x ~ so that £jl/£|| ~ \/k^{t±/t\\ 
result ~ £±A|| predicted by a 3D anisotropic NLaM.t3 



(2.17) 



Again this differs from the 



D. Effect of a parallel magnetic field 

We consider in this section the effect of an external magnetic field H = (0, H, 0) parallel 
to the planes. In second quantized form, the interplane hopping Hamiltonian is 

H± = -f± £ / dhe ie ^ )A{syds ^(v,l)i;(T,l') 
{i,i') J 

= -t± £ i/>\k\\ + (!-l')G,l')il>{k h l) (2-18) 

<i,j'>,k|| 



in the gauge A(0,0,-Hx). G = (G, 0, 0) with G = -eHd. In the second line of ( gig ) 



we have used a mixed representation by taking the Fourier transform with respect to the 
intraplane coordinate r. The effective action S[Q] of the NLaM can be simply obtained by 
calculating the particle-hole bubble II(q, uj v ) in the semiclassical (diffusive) approximation. 
To lowest order int±, we have IT = n^ + nW + n^ where IT^ ~ 27rJV 2 (0)r(l-|w I/ |r-Drgj[) 
(for \u v \t, Drq 2 <C 1) is the 2D result. Here oj v is a bosonic Matsubara frequency, llW and 

n(2) 

are given by (see Fig. H) 

n(1) = TT E G 2 {k h u n )G(k\\-8G,u n )G{k h u n+v ) 

k||,5=±l 

t 2 T 3 

= -8tcN 2 {0) 



l +w 2 r 2)l/2 ' 

n(2) = | E e^ w G(k| h o; n )6?(k|| - 5G,o; 7l )G(k| | ,a; n+ ,)G'(k|| - 5G,u; n+1 ,) 

k||,<5=±l 

= 87rN a (0) (1+ ^ )1/a cos( gi .d) , (2.19) 

for u) n oj n+v < and q\\,u u — > 0. L 2 is the area of the planes and 

G(k||,u n ) = {iu n + (i/2r)sgn(cu n ) - kj/2m + / u)" 1 (2.20) 

is the 2D one-particle Green's function (/i is the Fermi energy). We have introduced the 
characteristic magnetic energy uj c = VpG. The diffusion modes are determined by 

1 - (2vriV 2 (0)r)- 1 n(q, u v ) = \u u \r + Drq 2 + ^f 2 , ll2 sm 2 (q ± d/2) , (2.21) 

which yields the following interplane coupling in the NLaM: 
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S±[Q) = ~N 2 (0) — £ /d 2 r £ «QSK*)«Qfan(*) ■ (2.22) 

\ ' c J (1,1') inmafi 

Notice that we have retained only the diagonal part of the field Q since the magnetic field 
suppresses the interplane diffusion modes in the particle-particle channel. This is strictly 
correct only above a characteristic field corresponding to a "complete" breakdown of time 
reversal symmetry (this point is further discussed below). Eq. Q2.22p shows that the magnetic 
field not only breaks down time reversal symmetry but also reduces the amplitude of the 
interplane hopping. This quantum effect (important only when uj c t ^> 1) can be understood 
from the consideration of the semiclassical electronic orbits." Notice that such an effect 
cannot be described in the semiclassical phase integral (or eikonal) approximation for the 
magnetic field. The action S2D of the independent planes is not modified by the magnetic 
field since the latter is parallel to the planes. 

We can now apply the same RG procedure as in the preceding section. The initial (2D) 
stage of the renormalization is not modified by the parallel magnetic field, and the 2D/3D 
dimensional crossover is determined by 

(2 - 23) 

The crossover length L x (tj_,u> c ) therefore satisfies the scaling law 

ti 




_Ji ) 

1 + U^T 2 ) 1 /* J ' 



^2.24) 



where L x (t±) is the zero field crossover length obtained in the preceding section. 

In the 3D regime, the diffusion modes are 3D. The smallest volume corresponding to a 
diffusive motion is of order L 2 x d. This corresponds to a magnetic flux HL x d for a parallel 
field. Time reversal symmetry is "completely" broken when this magnetic flux is at least of 
the order of the flux quantum, i.e. when L x •> 1/G. It is easy to verify that this defines a 
characteristic field H such that oj c t <C 1. Above H , the diffusion modes are completely 
suppressed in the particle-particle channel. The effective action is then given by (|2.9| ) where 
only the diagonal part (i = j) of the field Q should be considered. D x and A x depend on H 
according to fl2.23j ). 

The critical value of the interplane coupling for H > H is obtained from g x = g' c where 
g' c is the critical value of the dimensionless conductance in the unitary case (no time-reversal 
symmetry). Thus, we have 

t { l\H)=t { f{l+u J y) l '\ (2.25) 

where ~ \g'J (TiV" 2 (0)^| D )] 1 / 2 < is the critical coupling for u c = in the unitary case. 
Close to the MIT, we have L x ~ £ 2 d (since g' c ~ 1) so that H is defined by u c ~ vpf&D, 
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which corresponds to an exponentially small value of the field: u c ~ T -i e -ak F i ^he ph ase 
diagram in the (t± — uo c ) plane is shown in Fig. |j. For H < Hq, the curve is not correct 
and should reach at H = 0. Therefore, we obtain that a weak magnetic field favors the 
metallic phase (t^(H) < t± for < H <C H ) while a strong magnetic field favors the 
insulating phase (t^(H) > for H > H ). This agrees with the results of Ref. [5] obtained 
from scaling arguments based on the weak localization correction to the Drude-Boltzmann 
conduct ivity.0 

In the metallic phase, the anisotropy of the correlation lengths is given by 




fi \ 1 



In the insulating phase, the anisotropy of the localization lengths is given by 



(2.26) 



(2.27) 



III. AUXILIARY FIELD METHOD 

The aim of this section is to recover the results of the RG approach by means of a 
completely different method. In order to study how the "correlations" are able to propagate 
in the transverse direction, we will study an effective field theory which is generated by a 
Hubbard Stratonovitch transformation of the interplane coupling Sj_[(3]. This methods bears 
some obvious similarities with standard approaches in critical phenomena.^ It is particularly 
useful when the low dimensional problem (here the 2D Anderson localization) can be solved 
(at least approximately) by one or another method. This kind of approach has already 
been used for weakly coupled ID Ginzburg-Landau mo delH or weakly coupled Luttinger 
liquids.0 While these studies were done in the high temperature (or disordered) phase, we 
start in this paper from the "ordered" (i.e. metallic) phase where the low-energy modes are 
Goldstone modes. 

Introducing an auxiliary matrix field ( to decouple S±, we rewrite the partition function 

as 

Z = J VQe'^ 1 ' 1 ' I d2rTt ^ J llM r ^ J pQ e -S 2D [Q]+2^ i /d2rTr[C i (r)Q i (r)] _ (g-Q 

The field ( should have the same structure as the field Q and therefore satisfies the condition 
£+ = C T ( T C = (. J xi ii is the inverse of the matrix Jxi,v = J_l(Si,i>+i + <5/,/'-i) where 
J± = (7r/4)iV2(0)t5_r (the matrix Jx is diagonal in the indices a,n,i). 

We first determine the value of the auxiliary field in the saddle point approximation. 
Assuming a solution of the form ij(C SF )Tn m ( T ) — Co^i,j^a,^n,m, we obtain the saddle point 
equation 

Co = J±(q± = 0) (uQ?Z(r)) SP , (3-2) 



S 



where J±(q±) = 2Jj_cos(q±d) is the Fourier transform of Jja,v- The average (• • -)sp should 
be taken with the saddle point action 

SsAQ] = S2d[Q] - 2 W rf 2 rTr[C SP g,(r)] - (3.3) 

i J 

Note that the saddle point value £ acts as a finite external frequency. The mean value 
of the field Q is related to the density of states and is a non-singular quantity. We havei 
{aQfnn ( r ))sp = sgn(w„) which yields 

Co = J±(q± = 0)sgn(u; n ) = 2Jj_sgn(u; n ) . (3.4) 

This defines a characteristic frequency for the 2D/3D dimensional crossover 

u x = — ^— \Q=4t 2 ,T. (3.5) 

We now consider the fluctuations around the saddle point solution. As in the standard 
localization problem, the Sp(2iV) symmetry of the Lagrangian is spontaneously broken by 
the "frequency" to Sp(iV) x Sp(iV).! The term Tr[f2Q] breaks the symmetry of S 2 d[Q]- Via 
the coupling Tr[£Q] between the fields ( and Q, it also breaks the symmetry of the effective 
action S[(] of the field (. Our aim is now to obtain the effective action of the (diffusive) 
Goldstone modes associated with this spontaneous symmetry breaking. Following Ref. || 
we shift the field according to ( — > ( + ( SP — 7rAT 2 (0)fi/4 and expand the action to lowest 
order in Q and (. The partition function becomes 

x J DQe-^ra+aE./^IfiWO.W] ; (3. 6 ) 

where we have introduce the 2D action 

S 2D [Q] =S 2D [Q] + ^N 2 (0)Y / J d 2 rTi[QQ l ( r )]-2Y / J d 2 rTr[C SP Q*(r)] 

= |tf 2 (0)Z>£ J d 2 rTr[V ll Q l (r)] 2 -2Y: J d'rTr^Q^)] . (3.7) 

S 2 r> is the action of the decoupled planes (t± = 0) at the finite frequency u x . To proceed 
further, we note that 

J pg e -^ra+2E l /^K I (r)O l W] = Z 2D e w ^ , (3.8) 

where W^lis the generating functional of connected Green's functions calculated with the 
action S 2 dS$ Z 2D is the partition function corresponding to the action S 2 n- We have 



W[C) = 2^2 / rf 2 rTr[0(r)(Q,(r))5 2 J 
i J 

+2J2J d 2 ri d 2 r 2 Y: «C2 B (n)«^S B (n,r2)i*C£»(r ! .) + ••• (3-9) 



ijnmaf3 
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whereS 

ijRinm( r h r 2) = (jiQlmn( r l)ijQ?nm( r 2)) S 2n ' (3.10) 



Using the saddle point equation ( |3.2|) , we obtain to quadratic order in the field ( and lowest 
order in Q the effective action 



SIC] = E / d 2 nd 2 r 2 £ ^J^) [j^t ~ r a ) - ^(n, r,)^] ,,c£(r 2 ) 

ijnmaf) 

~N 2 (0) E / rfVTr[fiJli,0(r)] . (3.11) 



For u n u m < 0, ijR, _ m (ri, r 2 ) is the propagator of the Goldstone modes of the action S^-dIv]. 
Its Fourier transform is given by 

'^ (q '' ) ^ (q '' ) = ^(o)(^g + ^ ) - < 3 ' 12 > 

where D x is the exact 2D diffusion coefficient at the finite frequency u x . Notice that the 
finite frequency ui x gives a mass to the Goldstone modes. The preceding equation defines 
the crossover length L x = (D x /uj x ) 1 / 2 . Since R(q\\ = 0) = 1/2 J±(q± = 0), the fluctuations 
of ( around its saddle point value are massless for u n u m < 0. On the other hand, it is clear 
that the fluctuations are massive for u n uj m > 0. Having identified the Goldstone modes 
resulting from the spontaneous symmetry breaking, we now follow the conventional NLcrM 
approach.! We suppress the massive fluctuations imposing on the field ( the constraints 
C 2 = C Sp2 an d Tr£ = 0. Rescaling the field ( in order to have ( 2 = 1 and introducing the 
Fourier transformed fields C(q), we obtain 

S[(] = Jl E(-T V) " 2i?(q||))Tr[C(q)C(-q)] - f iV 2 (0)Tr[QC(q = 0)] . (3.13) 

In the 3D regime, q\\ < 1/L X and q± ^ l/d, we can expand J± 1 (q±) ~ 2i?(qy) in lowest 
order in q\\ and q± to obtain 

SIC] = lN 2 (0) Y,(D x q 2 + 2ti C / 2 rgi)Tr[C(q)C(-q)] - \ iV 2 (0)Tr[fiC(q = 0)] . (3.14) 

Going back to real space and taking the continuum limit in the z direction (which introduces 
a factor l/d), we eventually come to 

S[(] = ^N 3 (0) J ^[^TrfVuC] 2 + 2^ 2 rTr[V,C] 2 ] - |iV 3 (0) J rf 3 rTr[^C] • (3.15) 

The cut-offs are A x = L~ l in the longitudinal directions and l/d in the transverse direction. 
Eq. ( p . 1 5 ) is similar to the action ( j2.9[ ) we have obtained in the RG approach. The 



only difference is that the crossover length is not defined in the same way. In the RG 
approach, L x ~ D(L x )/u x is defined via the length dependent 2D diffusion coefficient while 
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the auxiliary field method involves the frequency dependent 2D diffusion coefficient. We 
approximate the latter by@ 

DM= D p , (3.16) 
where uj v is a bosonic Matsubara frequency. This yields 

The critical interplane coupling is determined by N 2 (0)D x = N 2 (0)D(uj x ) ~ 1 which 
leads again to (|2.4j ). The crossover length L x is shown in Fig. [I]. For t± > there is an 
agreement with the RG approach. The reason is that in the weak coupling limit (g ^ 1), 
D(L) and D(uj u ) approximately coincide.0 Deep in the insulating phase, this is not the 
case and the two approaches give different results: L x 3> £ 2 d i n the RG approach while 
L x ~ £ 2 _d i n the auxiliary field method. This disagreement should not be surprising since 
neither one of the two methods is exact. In the RG approach, the dimensional crossover 
is treated very crudely since all the effects due to t± are neglected in the first (2D) stage 
of the renormalization procedure. However, this method should give qualitatively correct 
results. In particular, we expect the result L x ^> £ 2 d for t± <C t± to be correct. In the 
auxiliary field method, the effective action (|3.15| ) was obtained by completely neglecting 
the massive modes (more precisely sending their mass to infinity). In principle, these latter 
should be integrated out, which would lead to a renormalization of the diffusion coefficients 
appearing in (|3.15|) . The comparison with the RG result suggests that this renormalization 
is important in the insulating phase. 

The generalization of the preceding results in order to include a parallel magnetic field 
is straightforward. The interplane coupling is now given by ( |2.22| ). The auxiliary field ( has 
therefore only a diagonal part {i = j). The saddle point equation yields 

«(C SP )rl(r) = 5 a , p 6 n , m sgn(co n )^N 2 (0) + % 2)1/2 ■ (3-18) 
This defines the crossover frequency 

At 2 , t , 

^ H + ^)V < 3 - 19 > 

In the 3D regime, the massless fluctuations around the saddle point value are described by 
the action 

7T f r 2t 2 d 2 T i IT r 

S[(] = ~N 3 (0) j ^[^TrfVuC] 2 + (1 + ^ 2r2)1/2 Tr[V,C] 2 ] - ^(0) J d 3 rTr[QC] , (3.20) 

with the usual constraints on the field (. Here D x is the coefficient diffusion calculated at 
the magnetic field dependent frequency u x . Again we recover the result of the RG approach, 
the only difference coming from the definitions of L x and A x . 
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IV. CONCLUSION 



Using two different methods, we have studied the Anderson MIT in quasi-2D systems. We 
have found that the critical value of the single particle interplane coupling is given by fl2.4| ) . 
Apart from the factor 1/y/kpl, this result agrees with the diagrammatic self-consistent theory 
of Anderson localization and with estimates based on the weak localization correction to the 
Drude-Boltzmann conductivity. Nevertheless, it differs from recent numerical calculations 
according to which ~ In the weak disorder limit (kpl ^> 1), this latter result 

seems to us in contradiction with the scaling theory of Anderson localization since the 2D 
localization length £ 2 d is exponentially large with respect to the mean free path I. Because 
of the latter property, we indeed expect an exponentially small value of the critical coupling 

with respect to the elastic scattering rate 1 / r.lll For a very large 2D localization length 
^2D, it seems unlikely that the dimensional crossover and the MIT can be studied with 
numerical calculations on finite systems. The numerical calculations of Ref. |5] are done in a 
strong disorder regime (kpl ~ 1). In this regime, the exponential dependence of t± on kpl 
may be easily overlooked. 

We have also studied the anisotropy the correlation (localization) lengths in the metallic 
(insulating) phase and shown that it differs from the result predicted by a 3D anisotropic 
NLoM. The phase diagram in presence of a magnetic field was also derived: our approach 
formalizes and extends previous results obtained for weakly coupled chains.! 
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FIG. 1. Crossover length L x vs interplane coupling Solid line: RG approach, dashed line: 
auxiliary field method. We have used fe^Z = 50 and &d/1 = 5 10 4 (I and were considered as 
independent parameters). 
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FIG. 2. Anisotropy of the correlation (localization) lengths for t± > t± (t± < t±')- Solid line: 
RG approach, dashed line: auxiliary field method. The parameters are the same as in Fig. | with 
kpd ~ 1. 
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FIG. 3. Diagrammatic representation of II^ 1 ) (the symmetric diagram is not shown) (a) and 
n^ 2 ) (b). The solid lines represent the 2D Green's functions while the dashed lines denote single 
particle interplane hopping. 
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FIG. 4. Phase diagram in a parallel magnetic field. The position of (> ) is arbitrary. 
At very low field, uj c t ~ If&D, the curve is not correct since one should have t^(H = 0) = 
(see text). 
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